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We investigate the normal selfenergy and the mass enhancement factor in the Hubbard 
model on the two-dimensional square lattice. Our purpose in this paper is to evaluate the 
mass enhancement factor more quantitatively than the conventional third order perturbation 
theory. We calculate it by expanding perturbatively up to the fourth order with respect to 
the on-site repulsion U. We consider the cases that the system is near the half-filling, which 
are similar situations to high- Tc cuprates. As results of the calculations, we obtain the large 
mass enhancement on the Fermi surface by introducing the fourth order terms. This is mainly 
originated from the fourth order particle-hole and particle-particle diagrams. Although the 
other fourth order terms have effect of reducing the effective mass, this effect does not cancel 
out the former mass enhancement completely and there remains still a large mass enhancement 
effect. In addition, we find that the mass enhancement factor becomes large with increasing 
the on-site repulsion U and the density of state (DOS) at the Fermi energy p(0). According 
to many current reseaches, such large U and p(0) enhance the effective interaction between 
quasiparticles, therefore the superconducting transition temperature Tc increases. On the other 
hand, the large mass enhancement leads the reduction of the energy scale of quasiparticles, as 
a result, Tc is reduced. When we discuss Tc, we have to estimate these two competitive effects. 
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1. Introduction 

In recent decades, superconductivity in strongly corre- 
lated electron systems has been intensively investigated, 
motivated by the discoveries of various superconduc- 
tors including high- Tc cuprates,^ other transition metal 
oxides such as Sr2Ru04,^ heavy Fermion compounds'^ 
and molecular conductors.^ In order to describe the sys- 
tems theoretically, the Hubbard model has been often 
used as a minimal model including the on-site repulsive 
term. From the viewpoint of weak coupling, some au- 
thors calculated the momentum dependence of the effec- 
tive interaction using the fluctuation-exchange approxi- 
mation^^ and the perturbation theory,^ and then derived 
the anisotropic pairing state from the on-site repulsive 
interaction. In such studies, to obtain the realistically 
observable transition temperature Tc, one increases the 
magnitude of the bare interaction U. Large U enhances 
the effective interaction of quasiparticle, therefore Tc in- 
creases. However, it also enhances the effective mass,'' 
which leads the reduction of the energy scale of quasi- 
particles, as a result, Tc is reduced. When we discuss Tc 
quantitatively, we have to estimate these two competitive 
effects. 

In third order perturbation theory (TOPT),^ the third 
order terms for selfenergy cancel out completely for the 
case with the particle-hole symmetry, which occurs in the 
half-filled case of t' = 0. In this case, only one diagram in 
the second order contributes to the mass enhancement. 
For the doped case, the third order terms give the nega- 
tive contribution with large U, that is, the effective mass 
can not become large. To obtain more reliable results, 
we need to include higher order terms. The fourth order 



terms for the normal selfenergy have never been stud- 
ied before in lattice systems. They are composed of two 
kinds of the particle- hole diagrams, the particle-particle 
diagram and the vertex corrections. These particle-hole 
and particle-particle diagrams give the large mass en- 
hancement. The particle-particle diagram has the oppo- 
site sign to that in third order. As a result, when we 
calculated Tc in TOPT, we might have under-estimated 
the mass enhancement i.e. the reduction of the energy 
scale. 

For pairing interaction, Nomura et al. calculated the 
fourth order terms. ^ They showed the convergence of the 
pairing interaction is good for moderately strong U , when 
the system is near the half-filling. But they did not in- 
clude the normal selfenergy. In this paper, we calculate 
the mass enhancement factor by expanding the normal 
selfenergy up to the fourth order in the Hubbard model 
on the two-dimensional square lattice. We consider the 
cases that the system is near the half-filling. We hope 
that our quantitative investigation makes clear the es- 
sential mechanism determing Tc and explains the dif- 
ferences of observed Tc between the high- Tc systems^, 
YBa2Cu307_5 and La2-xSr2;Cu04 and so on. 

2. Formulation 

We consider the following single band Hubbard model 
on the two-dimensional square lattice. 



H — Ho + Hint , 



(1) 

(2) 
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(3) 

From the tight binding approximation, energy dispersion 
is given by 



^(fe) = — 2t ^cos kx + cos ky^ — At' cos kx cos ky 



(4) 



Here, t and t' are the nearest- and the next-nearest- 
neighbor hopping integrals respectively, and U is the 
on-site repulsion. Hereafter, we fix f = 1, and calculate 
physical quantities by changing the other variables t' ,U, 
electron filling n and the temperature of the system T. 
The bare Green's function is given as 



1 



(5) 



Here, fc is a shorthand notation as fc = (fc, ia;„) and a;„ = 
(2n+l)7rT (n : integer) is a fermion Matsubara frequency. 
The chemical potential /x is determined so as to satisfy 
the following equation. 



N 



(6) 



where; N is the number of the lattice sites. The bare sus- 
ceptibility Xo(<?) is given as 

k 

where g is a shorthand notation as q = {q,i^m) and 
Oto = ^rmrT (m : integer) is a boson Matsubara fre- 
quency. 

The normal selfenergy is expanded perturbatively in 
U as follows. 



S(fc) = S(2)(fc) + S(3)(fc) + S(4)(fc) + • • • 



(8) 



here S(")(fc) is simply proportional to U" because there 

is no momentum dependence in the on-site repulsion U. 
We evaluate the normal selfenergy up to the fourth or- 
der with respect to U. The diagrammatic and analytic 
expressions are given in Appendix A. 

In this paper, we mainly discuss the mass enhancement 
factor z~^{k), which is given as 

9SR(fe,w) 



-i(fc) = l 



duj 



(9) 



Here, SR(fe,w) is real frequency retarded selfenergy, 
which is obtaind by performing the analytic continua- 
tion of S(fc, iun) from the upper half plane with use of 
Pade approximation. 

We take 64 x 64 fe-meshs for the first Brillouin zone and 
512 Matsubara frequencies in the numerical calculations. 
We fix the temperature of the system T = 0.01. The 
other details of the calculation are given in Appendix A. 




Fig. 1. The Fermi surfaces in the case of t = 1, t' = —0.15 and 
T = 0.01. The solid (dashed) line is the result for n = 0.9 (n = 
1.1). The symbols F, D and X represent the symmetry points 
(0,0), (0,7r) and (tt, tt) respectively. 



3. Results 

We consider the cases that the system is near the half- 
filling. These are similar situations to high-Tg cuprates. 
We fixed the parameters t = l,t' = —0.15 and T = 0.01. 
In Fig.l, we show the Fermi surfaces for n = 0.9 and 
n = 1.1. As we set the half-filling n = 1, the elec- 
tron filling n = 0.9 (n = 1.1) represents the hole-doped 
(electron-doped) case. We can sec the (tt, tt) nesting, 
which yeilds the strong antiferromagnetic (AF) fiuctu- 
ations. The band width of these dispersions equals 8.0. 
Figure 2 shows the density of state in these two cases. 
The DOS at the Fermi level for n = 0.9 is larger, bacause 
when we introduce the negative t' at the half-filling, the 
level of the van Hove singularity falls down. We can see 
the strong peak of AF fiuctuations at (7r,7r) in Fig. 3. 
The peak is stronger in the hole-doped case because of 
the strong nesting property and the large p{0). 

We have calculated the mass enhancement factor 
z~^{k) in those two cases n = 0.9 and n = 1.1. First, we 
show the result for the hole-doped case n = 0.9 in Fig. 4. 
By introducing the fourth order terms, we obtain the 
large efii'ective mass at the Fermi surface, which is almost 
twice as large as that calculated in TOPT. This result 
solves the problem of the reduction of z^^{k) due to the 
third order terms as we mentioned before. We can define 
the value as the average of z^^(A;) at the Fermi point 
near (0,7r) and that near (7r/2,7r/2). In Fig. 5, we show 
the J7-dependence of z~^ in the case of n = 0.9. With in- 
creasing U, the mass enhancement factor increases. From 
this calculation, we can consider to be a function of 
U and estimate the coefficients of f/", such as 

^-1 = 1-1- 0.0564C/^ - 0.0048[/^ -|- 0.0063i7^. (10) 

Next, we show the result for n = 1.1, in Fig. 6. We can 

also see the large mass enhancement, but it is relatively 
low. This result is related to the fact that p{0) is lower 
than that in the case of n = 0.9. In Fig. 7, we show the 
[/-dependence of and estimate z~^ to be 

z-'^ = 1 + 0.0273t/^ - O.OOllt/^ + 0.001417'' (11) 
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Fig. 2. The density of state in the case of t = I, t' = —0.15 
and T = 0.01. The dotted line and the zero energy represent the 
Fermi level for n = 0.9. And the dashed line represents the Fermi 
level for n = 1.1. 
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Fig. 4. The mass enhancement factor 2~^(fc) in the case of n = 
0.9 and 1/ = 3. The solid, dotted and dashed lines are the results 
in the calculations up to the second, third and fourth order, 
respectively. 2~^(fc) = 1 means the mass of the free electrons. 



Xo(q,0) 

0.8 




Fig. 3. The bare spin susceptibilities Xo(Qi 0) in the case of t = 1, 
t' = -0.15 and T = 0.01. The solid (dashed) line is the result 
for n = 0.9 (n = 1.1). 




Fig. 5. The (7-dependence of the mass enhancement factor 
in the case of n = 0.9. The solid, dotted and dashed lines are 
the results in the calculations up to the second, third and fourth 
order, respectively. is the average of z~^(fc) at the Fermi 

point near (0,7r) and that near (7r/2, 7r/2). 



When we set a large U value, we can also obtain the large 
mass enhancement. 

As we discuss in Appendix B, roughly speaking, the 
contribution from all the fourth order terms is about half 
as much as the one from the diagram (4A) as the result 
of the cancellation in the hole-doped case (n = 0.9). In 
the electron-dope case (n = 1.1), the cancellation is more 
remarkable, and the ratio of the coefficient of to that 
of in eq.(ll) is smaller than that in cq.(lO). 

4. Discussion 

In this paper, we have investigated the the mass en- 
hancement factor perturbatively up to the fourth order 



in the Hubbard model. We obtained the large mass en- 
hancement on the Fermi surface compared with the re- 
sult in TOPT, and it becomes large with increasing I] 
and p(0). 

In this section, we discuss the relation between mass 
renormalization and T^. Some anthers have described the 
Fermi liquid theory for the strongly correlated electron 
systems ''.^^ Based on this idea, electrons form the quasi- 
particles with effective mass and finite lifetime due to 
the effect of the interaction U. Therefore, we have to 
start with the estimation of the renormalization factor 
z, then calculate the momentum dependence of the effec- 
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Fig. 6. The mass enhancement factor z~^{k) in the case of n = 
1.1 and (7 = 3. The soUd, dotted and dashed lines are the results 
in the calculations up to the second, third and fourth order, 
respectively. z~^{k) = 1 means the mass of the free electrons. 
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Fig. 7. The (7-dependence of the mass enhancement factor 
in the case of n = 1.1. The solid, dotted and dashed lines are 
the results in the calculations up to the second, third and fourth 
order, respectively. 

tive interaction between the renormalized quasiparticles. 
In that scheme, the renormahzed parameters are given 
as t^''^ = zt, t'^^^ = zt', = zU, where the index 

(r) denotes renormahzation due to the normal selfenergy 
correction and z is the renormahzation factor. Therefore, 
we consider that the transition temperature in the rcnor- 

(r) 

maUzed scheme Tc is given as 

T^"'' = zT, (12) 

This relation show the reduction of Tc due to small z 
as we mentioned before. Taking our calculation into ac- 
count, this tendency is considered to be more marked by 



introducing the fourth order terms. We have to perform 
the calculation of Tc based on the fourth order selfenergy 
correction. This is further investigation in the future. 

Here we would like to stress the following point. The 
anisotropic part of quasiparticle interaction determines 
the superconductivity. On the other hand, the mass en- 
hancement factor is determined by total interaction in- 
cluding isotropic interaction. Therefore, in principle, we 
have independent two parameters which determine the 
transition temperature. The two parameters are deter- 
mined by details of electronic structures, depending on 
the Fermi surface and d or / electron numbers. 

5. Conclusions 

We investigate the mass enhancement factor by ex- 
panding the normal selfenergy perturbatively up to the 
fourth order with respect to U in the Hubbard model. 
We consider the cases that the system is near the half- 
filling on the two-dimensional square lattice. These are 
similar situations to high- Tc cuprates. As results of the 
calculations, by introducing the fourth order terms, we 
obtain the large mass enhancement on the Fermi surface 
compared with the result in TOFT. This is mainly due to 
the fourth order particle-hole and particle-particle terms. 
Although the other forth order terms have effect of re- 
ducing the effective mass, this effect does not cancel out 
the former mass enhancement completely and there re- 
mains still a large mass enhancement effect. In addition, 
wc find that the mass enhancement factor becomes large 
with increasing the on-site repulsion U and the density 
of state at the Fermi energy p(0). According to many 
current reseaches, such large U and p{0) enhance the ef- 
fective interaction between quasiparticles, therefore Tc 
increases. On the other hand, the large mass enhance- 
ment leads the reduction of the energy scale of quasipar- 
ticles, as a result, Tc is reduced. When we discuss Tc, 
we have to estimate these two competitive effects. This 
is the important unified theory for all the strongly cor- 
related electron systems from cuprates to heavy fermion 
compounds. 
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Appendix A: Details of Calculation 

The normal selfenergy is expanded perturbatively in 
U as follows, 

E(fc) = (fc) -I- e(3) (fc) + (fc) + ■ • • , (A-1) 

here I]''"-'(fc) is simply proportional to t/". These pertur- 
bation terms are represented by using the bare Green's 
function G^''-'(fc), and the bare susceptibilities, 

fe 

'^"('?) = -§EG^"^(9"fc)G(°)(fc). (A.3) 
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Fig. A-1. The Feynman diagrams of the normal selfenergy up to 
the fourth order. Solid and dashed lines correspond to the bare 
Green's function and the interaction, respectively. 



The diagrams for the normal selfenergy up to the fourth 
order are shown in Fig.A-1, and the analytic representa- 
tions for them are given in the followings, where k is the 
shorthand notation of the external momentum and fre- 
quency, and Pi's (ft's) are those of the internal momenta 
and frequencies for fermions (bosons). 

The second order term 

S(2)(fc) = C/2^^ Xo(g)GW(fc-g). (A.4) 
The third order term 



S = f,3i: ^ ^^(^)2 ^(0)(fc _ (A.6) 



S(3vc)(fc) = u'L Yl Mlf - k). (A.7) 



N 



The fourth order term 



(fc) = 5](4A) + (fc) + . . . + S(4L) (A;), (A-8) 

where 



S(^A)(fc) = 




Xo{qf G(o)(fc-g), 


(A-9) 






Xo(9)' G(°nfc-«), 


(A-IO) 


E(4C)(A:) = 


N ^ 

9 


Mlf G(°)(g-fc), 


(A-11) 




p 







E(4H)(A:) 



(A-12) 



(A-13) 



E(4F)(fc) = C/^^^ ^o(/c+p) G(°)(p)2 E(2)(p), 



(A.14) 



s(^^Hfc) = t^^(^)'EExo(?i)xo(g2) 



91 92 



X G(°)(fc-gi) G(°)(fc-g2) 

xGW(fc- 91-92), (A-15) 



91 92 



xGW(fc-9i) GW(g2-A;) 

xGW(gi + g2-fc), (A-16) 

E(«)(fc)= - f/^(|)^^^xo(gi) MQ2) 

91 92 

xG(°Hfc-<Zi) G(°)(g2-fc) 
xG(°)(gi+g2-fc), (A-17) 



s(")(fc) = t/^(^)^EEExo(3-Pi-P2) 

9 pi P2 

xG(")(pi) G(")(p2) G(")(g-pi) 
xG(")(g-p2) G^^H?-^), (A-18) 
E(4«)(fc) = [/4(|;)^EEE 'Ao(9+Pi+P2) 



9 Pi P2 
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xG^°\q+P2) G^'^Hk-q), (A- 19) 

q Pi P2 

xG(°)(g+p2) G(")(A,--g). (A-20) 

We take 64 x 64 fe-meshs for the first Brillouin zone 
and 512 Matsubara frecencies in the numerical calcula- 
tions. Wc evaluate the second order term, the third order 
terms, the fourth order (4A), (4B), • • • , (4E) and (4F) 
terms by using the fast Fourier transformation (FFT) 
algorithm because the summations included in them all 
have convolution forms. In the cases of the fourth order 
(4G), (4H) and (41) terms, we fix the external momentum 
and frequency k. Then, the full summations with respect 
to the internal momenta and frequencies are performed 
by using the FFT algorithm. 

The evaluation of the fourth order (4J), (4K) and (4L) 
diagrams is very difficult. The summations of momenta 
and frequencies for them are not reduced to convolution 
forms even if wc fix the external momenta k, so the es- 
timated computation time is not realistic. Therefore, we 
take the following approximation. First, we notice that 
these diagrams have such forms as only one bare Green's 
function contains the enternal momenta k, then we sep- 
arate it and write the rest as A{q). For example, in the 
case of (4J), 

S(«)(A;) = C/'^E^W G^^Hq-k) (A.21) 
q 

Ml) = {j^YJ2Yl Xo{q-Vi -P2) 

Pl P2 

xG(°)(pi) G(°)(p2) G(°)(9-pi) 
xG^°\q-p2). (A-22) 

The equation (A-21) has convolution form, therefore the 
remaining problem is the calculation of A{q). We per- 
form the Fourier transformation and consider the follow- 
ing quantity, 

Air)=J2 A{'l)e""'. {A-23) 

q 

This quantity A{r) is reduced to 
q' 

Pl 

^ § E (^^2) G(°) {q' + P2) e^'''^ (A.24) 

P2 

where r is the shorthand notation r = (r, r), and we 
have defined the variable q' = q — pi — p2- We can see 



Author Name 

that the summations in eq.(A-24) can be performed with 
use of the FFT algorithm for fixed r. We calculate A{r) 
for only 9x9 r-region around \r\ = 0, because with 
increasing |r|, A(r) converges rapidly. This is a similar 
way which Nomura et al. estimate the non-Parquet 'en- 
velope' diagram.^ This cut-off of the large \r\ components 
means that in q-space the small structure of A{q) is ne- 
glected. Here, we note that although the shape of A{q) 
becomes rather smooth, this approximation hardly affect 
the normal selfenergy quantitatively. As for calculating 
in r-space, we reduce the r-mesh to 256, then we fill the 
middle points with the averages. We have verified that 
this reduction also has no quantitative problem. Finally, 
we can obtain ^(''•'^(fc) by calculating the convolution of 
A{q) and G(°) {q - k) with use of the FFT algorithm. 

Appendix B: Contribution from Each Diagram 
in Fourth Order 

In this section, we briefly discuss the contribution from 
each diagram in fourth order. The RPA diagrams (4A) 
and (4B) give the large mass enhancement, and so is 
the particle-particle diagram (40). Only these three di- 
agrams enhance the effective mass, and then the other 
nine diagrams (4D), • • • , (4K) and (4L) reduce it. If we 
neglect the lattcrs, we obtain about six times larger co- 
efficient of in eq.(lO), that is too overestimating. The 
diagrams (4G), (4H) and (41) give large negative con- 
tribution to the mass enhancement, which mainly can- 
cel the formers' positive contribution however not com- 
pletely. On the other hand, the one from the diagrams 
(4D), (4E) and (4F) is very small. And the fact does not 
change even if we collect this type of diagram up to in- 
finite order. In the case of the diagram (4J), (4K) and 
(4L), their contribution is about half as much as that 
from the diagram (4G). This is considered to be related 
to the fact that there is one extra summation which can 
not be reduced to xo or 0o- As the result of the cancel- 
lation, roughly speaking, the contribution from all the 
fourth order terms is about half as much as the one from 
the diagram (4 A) in the hole-doped case (n = 0.9). In 
the electron-dope case (n = 1.1), the cancellation is more 
remarkable, and the ratio of the coefiicient of to that 
of U"^ in eq.(ll) is smaller than that in eq.(lO). 
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